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GEOMETRIC  MODEL  FOR  TRACKER-TARGET 
LOOK  ANGLES  AND  LINE  OF  SIGHT  DISTANCE 

By 

DANIEL  T.  LAIRD 


ABSTRACT 

To  determine  the  tracking  abilities  of  a  Telemetry  (TM)  antenna  control  unit  (ACU) 
requires  ''truth  data"  to  analyze  the  accuracy  of  measured,  or  observed  tracking  angles.  This 
requires  we  know  the  actual  angle,  i.e.,  that  we  know  where  the  target  is  above  the  earth.  The 
positional  truth  is  generated  from  target  time-space  position  information  (TSPI),  which  implicitly 
places  the  target’s  global  positioning  system  (GPS)  as  the  source  of  observational  accuracy.  In 
this  paper  we  present  a  model  to  generate  local  look-angles  (LA)  and  line-of-sight  (LoS)  distance 
with  respect  to  (w.r.t.)  target  global  GPS.  We  ignore  inertial  navigation  system  (INS)  data  in 
generating  relative  position  at  time  T;  thus  we  model  the  target  as  a  global  point  in  time  relative 
to  the  local  tracker’s  global  fixed  position  in  time.  This  is  the  first  of  three  companion  papers  on 
tracking  This  is  the  first  of  three  companion  papers  on  tracking  analyses  employing  Statistically 
Defensible  Test  &  Evaluation  (SDT&E)  methods. ** 


KEY  TERMS 

TM  antenna,  ACU,  Latitude,  Longitude,  Altitude,  Look-Angle,  Line-of-Sight,  TSPI, 
GPS,  Local  Plane,  Spherical  Trigonometry,  Law  of  Spherical  Cosines;  Law  of  Spherical  Signs. 


INTRODUCTION 

The  target  TSPI  file  is  our  truth  source.  It  contains  GPS,  and  often  INS  data.  Our  task  is 
to  translate  target  global  latitude  (TLAT),  longitude  (TLNG)  and  altitude  (TALT)  at  time  (T)  to 
local  tracking  ‘look-angles’  of  azimuth  (AZ),  elevation  (EL)  and  line-of-sight  distance  (LoS)  at 
time  T.  Time  correlates  target  position  (SLAT,  SLNG,  SALT)  with  a  fixed  location  as  tracker 
site  origin.  Below  we  develop  a  set  of  global-to-local  transformations  between  tracker  and  target 
at  time  T  using  derived  target  angles  and  target/tracker  angle  differences: 

=  TLAT  -  SLAT,  (la) 

AZ  =  TLNG  -  SLNG,  (Ib) 

AH  =  TALT  -  SALT.  (Ic) 

Our  geometric  model  is  merely  one  of  a  set  of  tools  for  we’ve  developed  for  estimating 
tracking  efficiency.  It  serves  as  our  primary  analysis  space*^.  The  purpose  of  our  geometric 
model  is  to  estimate  ACU  tracking  errors  and  mode  probabilities.  Error  is  not  exclusive  to  an 
ACU,  but  rather  extends  to  the  scope  of  the  test  scenario.  Target  inertial  data  is  useful  for 
modeling  tracking  angle  errors  due  to  tracker-target  antenna  alignment.  A  second  paper  focuses 


"  Rf.  1  &  2. 

Rf.  3. 
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on  modeling  tracking  error  which  employs  GPS  and  INS;  a  third  on  modeling  autotracking  mode 
of  an  ACU  employing  receiver  gain  (AGC)  and  antenna  scanning  controls. 

LOCAL  AZIMTUH  AND  ELEVATION  MODELS: 

The  flat  azimuth  model  is  shown  in  figure  1.  Basically,  two  local  areas  on  an  earth 
tangent  plane  (an  earth  local  frame),  with  each  locally  flat  plane  a  distance,  ‘d’  apart  -  the 
calculated  earth  distance,  which  is  d  =  R  x  cos"^(L),  where  L  denotes  the  earth  surface  arc 
between  the  two  points.  LoS  designates  the  distance  between  points  at  altitude.  In  the  distance 
equation,  the  earth  radius  is  approximated  as  R  6378137m,  and  0  is  calculated  via  projections 
of  latitude,  longitude  and  longitudinal  differences  via  the  spherical  geometry. 


Our  model  formulae  are  derived  from  spherical  law  of  cosines  (cf  below).  There  are 
several  estimates  of  arc  length,  ‘d’  on  sphere.  We  employ  two  (cf  figure  1.1  below): 


cos(d)  =  sin((t)i)  x  sin((t)2)  +  cos((t)2)  x  cos((t)i)  x  cos(AA,),  (la) 

L  =  Rxacos(d).  (lb) 

Substituting  d  ^  L/R  yields  a  more  accurate  close-distance  formula: 

2  2  2 

sin  (l/il)  =  haversin  (1)  =  sin  -I-  cos((t)i)  x  cos((t)2)  x  sin  (ViAX),  (2a) 

L  =  R  X  2asin(Vl).  (2b) 
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The  haversine  is  more  aeeurate  for  smaller  arc  distances  (angles);  where  smaller  is 
relative  w.r.t.  R.  Mapping  the  plane  of  figure  1  onto  a  sphere,  we  get  a  grid  shown  in  figure  1 . 1 
below.  In  this  figure,  from  a  perspective  above  the  plane  (j)  designates  latitude,  X  longitude  and  L 
designates  the  projection  of  LoS  onto  the  spherical  plane. 


Figure  1.1  Mapping  of  Plane  in  Figure  1  onto  a  Sphere 

The  azimuthal  calculations  are  based  on  the  law  of  cosines  and  sines  for  a  spherical 
coordinate  system: 

0  =  (ISO/tt)  X  atan2((p)  (3a) 

(p  =  (sin(A^)  X  cos((|)i),  sin({|)i)  x  cos((t)2)  -  cos({|)i)  x  sin({|)2)  x  cos(AA,)).  (3b) 


Figure  2  Planar  Geometric  Elevation  Model 

A  flat  or  planar  geometric  elevation  model,  shown  in  figure  2  is  based  on  the  difference 
in  altitudes  and  the  LoS  or  L: 
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\|/  =  (I8O/71)  X  asin(AH/LoS)  =  atan(AH/L). 


(4) 


The  arctangent  is  preferred  as  AH  >  LoS  is  possible  during  flight.  The  elevation  angle  is 
the  angle  of  altitude,  \|/.  The  LoS  and  azimuth  angle  formulae  are  derived  below.  First  we  prove 
the  theoretical  ‘truth  source’,  i.e.,  the  laws  of  spherical  trigonometry. 

Law  of  Spherical  Cosines  (LSC): 

The  Law  of  Spherical  Cosines  is  one  of  two  fundamental  formulae  used  as  our  navigation 
truth,  so  we  shall  prove  it  and  its  companion,  the  Law  of  Spherical  Sines.  Of  course,  this 
theoretical  ‘truth  source’  compares  with  observed  ‘truth  source’  that  is  GPS  TSPI.  Using  the 
triangle  in  figure  3  the  LSC,  or  the  so-called  cosine  rule  of  sides  is  represented  by  the  formula: 

cos(c)  =  cos(a)cos(b)  -  sin(a)sin(b)cos(C).  (5) 


Figure  3  Spherical  Reference  Points 


LSC  PROOF: 

This  proof  is  attributed  to  Romuald  Ireneus  (1997).  We  use  the  unit  sphere,  as  the  radius 
will  drop  out  of  the  calculations.  For  unit  vectors  u,  v,  w  from  center  of  the  sphere  to  the  comers 
of  the  triangles,  the  inner-product  of  one  vector  on  the  other  yields  the  cosine  angle  subtended  by 
arc  length:  =  1  x  0.  The  cosine  on  the  sphere  is  still  the  inner-product  of  the  subtending  units. 
cos(a)  =  u  •  v;  cos(b)  =  u  •  w;  cos(c)  =  w  •  v.  (6) 

We  define  tangents  perpendicular  to  each  unit: 

ta  =  V  -  u(u  •  v)/||  u  -  v(u  •  v)  II  =  (v  -  ucos(a))/sin(a),  (7a) 

tb  =  =  (w  -  ucos(b))/sin(b).  (7b) 

Therefore,  the  ‘law  of  arcs’  is: 
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cos(C)  =  ta  •  tb  =  (v  -  ucos(b))(w  -  ucos(b))/sin(a)sin(b)  (8a) 

=  (cos(c)  -  cos(a)cos(b))/sin(a)sin(b) 

cos(c)  =  cos(a)cos(b)  +  sin(a)sin(b)cos(C).  (8b) 

Similarly,  ‘the  law  of  angles’  is: 

cos(A)  =  -cos(B)cos(C)  +  sin(B)sin(C)cos(a).  (8c) 

QED. 


Law  of  Spherical  Sines  (LSS): 

Again  using  the  triangle  from  figure  3  above,  the  LSS  or  so-called  rule  of  sines  is 
represented  by  the  formula: 

sin(A)/sin(a)  =  cos(B)/cos(b)  =  sin(C)/sin(c).  (9) 


Figure  4  Spherical  Model  Frame 


LSS  PROOF: 

From  Pythagoras’  law  we  know  that  for  any  angle  (arc  length),  A 

sin^(A)  =  1  -  cos^(A).  (10a) 

Using  the  LSC  from  above: 

sin^(A)  =  1  -  (cos(a)  -  cos(b)cos(c))/sin(b)sin(c))^.  (10b) 


A  little  manipulation  yields: 

sin(A)/sin(a)  =  cos(a,b,c)/sin(a,b,c) 
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(10c) 


where 


2  2  2  !/ 

cos(a,b,c)  =  [1  -  cos  (a)  -  cos  (b)  -  cos  (c)  +  2cos(a)cos(b)cos(c))]  ^ 

sin(a,b,c)  =  sin(a)sin(b)sin(c). 


(lOd) 

(lOe) 


A  permutation  on  pairs  (A, a),  (B,b),  (C,c)  proves  the  rule.  This  proof  is  attributed  to  Isaac 
Todhunter  circa  1863.  We  will  use  these  spherical  trigonometric  laws  to  construct  the  azimuth 
formulae. 

AZIMUTH  AND  LOS  PROOFS: 

For  instructional  insight  we  derive  azimuth  in  two  ways:  the  first  using  vectors,  the 
second  via  the  spherical  trigonometry  laws  derived  above.  Our  analyses  calculations  employ  the 
trigonometric  formulae,  without  taking  account  of  the  earth’s  elliptical  eccentricity  -  our 
relatively  small  distance  w.r.t.  R  between  target  and  tracker  permit  this.  The  spherical  model 
frame  is  shown  in  figure  4.  The  error  due  to  ignoring  eccentricity  is  insignificant  for  the 
precision  of  our  calculations.  To  derive  azimuth  via  vectors  we  create  unit  vectors  on  the  sphere 
referenced  to  a  cartesian  coordinate  system,  e.g.: 

ek  =  (cos((|)k)cos(?Lk),  sin((|)k)cos(;Lk),  sin((|)k))  <»  (ei(k),  e2(k),  e3(k)).  (1  la) 

I  found  this  approach  in  an  earthquake  seismology  class  note®,  and  believe  the  derivation  and  an 
alternative  enlightening,  so  include  both  here  for  completeness.  For  two  coordinates  we  have  an 
arc  length  (on  unit  sphere)  given  by  the  inner-product: 

COS(A)  =  Cl  •  62.  (11b) 

Thus  the  arccosine  will  provide  angle.  Simple!  Although  mathematically  simple,  when  angles 
are  small,  round-off  errors  emerge  and  grow  as  cos(A)  — >  1.  So,  we  derive  an  equivalent 
alternative  formula  that  employs  the  sine  and  cosine  to  create  the  tangent;  and  though  more 
complicated,  is  numerically  more  accurate  for  computer  computation.  Figure  5  shows  a  way  to 
construct  the  sine  and  cosine  of  the  bisected  angle: 


ei  -  62 


Figure  5a  Sine  and  Cosine  of  Bisected  Angle 


We  derive  the  so-called  haversine  representation  of  LoS  from  the  vectors  of  figure  5a: 
sin(V2A)  =  Vij  61  -  62I,  cos(1/2A)  =  Vij  61  +  62I  —>■  sin(A)  =  2sin(V2A)cos(V2A).  (11c) 


So, 


A  =  atan2(sin(A),cos(A)), 


(lid) 


Rf.  4  &  8. 
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(lie) 


Where,  from  11b  and  the  sine-cosine  relations  we  have: 
cos(A)  =  ei  •  62,  sin(A)  =  Vij  ei  -  62]  x  |  ei  +  62!. 


N 


Figure  5b  Grid  for  Spherical  Trigonometry  Laws 
Now  we  make  the  simple  derivation  using  the  spherical  trigonometry  laws.  Using  the 


grid  of  figure  5b‘*  yields: 

sin(0)/sin(!/27i  -  {|)i)  =  sin(AA,)/sin(A)  =  sin(27i  -  0*)/sin(!/27i  -  (|)2),  (12a) 

cos(7i  -  (|)i)  =  cos(A)cos(7r  -  {|)2)  +  sin(Ai2)sin(7r  -  {|)2)cos(0).  (12b) 

These  are  rewritten  as 

cos((|)2)sin(A)sin(0)  =  cos((|)i)cos((|)2)sin(A^),  (12c) 

cos({|)2)sin(A)cos(0)  =  sin((|)i)  -  sin((|)2)cos(A).  (12d) 

Therefore, 

0  =  (180/7r)  X  atan2(sin(0),cos(0))  =  atan2((p),  (13a) 

cp  =  cos({|)i)cos((|)2)sin(A^)/(sin((|)i)  -  sin((|)2)cos(A^)).  (13b) 


Azimuth  is  the  angle  of  (13b)  which  is  equivalent  to  (14a)  below.  Reversing  the 
arctangent  arguments  supplements  perspective,  i.e.,  we  move  the  look-angles  from  tracker-to- 
target  to  target-to-tracker. 

QED. 


For  the  haversine  note  that 

haversin(0)  =  sin^(0/2);  cos(0)  =  1  -  2sin^(0/2)  =  1  -  2haversin(0).  (14a) 


^  Ibid. 
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Now  substitute  this  form  of  cos(0)  and  employ  the  identity: 
eos(a  -  b)  =  eos(a)cos(b)  +  sin(a)sin(b), 


(14b) 


into  equation  (8b)  above  to  obtain  the  law  of  haversines: 

h  =  haversin(L/R)  =  haversin(A(t))  +  sin(a)sin(b)eos(A>.),  (14e) 

L  =  2Rasin(\h) .  ( 1 4d) 

Where  ‘AAcj)  =  a  -  b  and  VilsX  =  C,  the  latitude  and  longitude  differenees  of  figure  4. 
QED. 


Observed  Elevation  Angle(ACU) 


Calculated  Elevation  Angle(TSPI) 


in 
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Time[seconds) 


Observed  Azimuth  Angle(ACU) 


Calculated  Azimuth  Angle(TSPI) 


Tim0(secQnds) 


Tim0(seconds) 


Figure  6a  Observed  and  Modeled  Angles 
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Observed  Elevation  Angle(ACU} 


Observed  Azimuth  Angle{ACU} 


Figures  6a  and  6b  show  observed  (track  measured)  and  theoretical  (calculated  from  TSPI) 
azimuth  and  elevation  angles.  The  sample  times  are  seconds  from  test  start  time,  T  =  To  =  0.  It 
is  clear  that  between  6k  to  8k  seconds  that  tracking  error  increased,  significantly  w.r.t.  to  the  rest 
of  the  profde.  From  the  whole  angle  tracking  profde  we  conclude  that  overall  the  measured  was 
quite  accurate  w.r.t.  the  calculated,  or  theoretical  tracking  profile.  Figure  6b  shows  the  error 
superimposed  on  the  angle  means,  which  is  a  numerical  verification  of  our  conclusion.  A  sample 
of  observed  and  calculated  angles  is  shown  in  figure  7. 
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^  Az.EI.Angles.sample.txt  -  Notepad 


File 

Edit  Format 

View  Help 

>  AZ 

El. Angles [1000:1010,] 

> 

Ti  me 

AZ. Actual 

El 

Actual 

AZ.  Cal  c 

El  .Calc 

AZ. Err 

El . Err 

1000 

13: 50:Z7 

40. 83Z 

Z.  819 

37. 3Z431 

3.131660 

-3. 507695 

0. 31Z660 

1001 

13: 50:Z8 

40. 6Z3 

Z.  835 

37.1Z1Z6 

3.14115Z 

-3. 501738 

0. 30615Z 

lOOZ 

13: 50:Z9 

40. 398 

Z.  846 

36. 95551 

3.154514 

-3.44Z488 

0. 308514 

1003 

13:50:30 

40. Z16 

Z.  846 

36. 77Z59 

3.169780 

-3.44341Z 

0. 3Z3780 

1004 

13:50:31 

40. 008 

Z.  846 

36. 57514 

3.18Z340 

-3.43Z86Z 

0. 336340 

1005 

13:50:3Z 

39.788 

Z.85Z 

36.41386 

3.195Z45 

-3.374143 

0.343Z45 

1006 

13:50:33 

39. 579 

Z.  879 

36. Z1960 

3. Z07480 

-3. 359401 

0. 3Z8480 

1007 

13:50: 34 

39. 376 

Z.  879 

36. 04359 

3. Z15Z07 

-3. 33Z410 

0. 336Z07 

1008 

13:50:35 

39.189 

Z.  879 

35. 8871Z 

3. ZZ7684 

-3. 301878 

0. 348684 

1009 

13:50:36 

38. 991 

Z.  879 

35. 69793 

3. Z36094 

-3. Z9307Z 

0. 357094 

1010 

13:50:37 

38. 799 

Z.  91Z 

35.  54035 

3. Z48Z58 

-3.  Z:58650 

0. 336Z58 

>  AZ. 

, El. Angles [2000: 2010,] 

Ti  me 

AZ. Actual 

El, 

.Actual 

AZ.  cal  c 

El .  cal  c 

AZ. Err 

El . Err 

zooo 

14 : 07 : 07 

6.768 

3.483 

5. 813365 

3. 597Z84 

-0.954635 

0.114Z84 

ZOOl 

14:07:08 

6.  741 

3.47Z 

5. 797938 

3. 59547Z 

-0. 94306Z 

0.1Z347Z 

zooz 

14:07:09 

6.  730 

3.47Z 

5. 779Z60 

3. 593769 

-0. 950740 

0.1Z1769 

Z003 

14:07:10 

6.  708 

3.47Z 

5. 763937 

3. 591998 

-0. 944063 

0.119998 

Z004 

14:07:11 
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Figure  7  Observed  and  Calculated  Angles 


SUMMARY 

To  determine  tracking  accuracy  we  need  both  a  measured  data  and  theoretical  truth 
source  to  compare  observed  angles  with  expected  angles.  The  equations  used  to  compare 
measured  with  expected  angles  are: 


Line  of  Site  Distance: 

LoS  =  2R  X  asin(A/h),  R  ~  6378137m 

(2a) 

h  =  haversin(A(|))  +  sin((|)i)sin((|)i)cos(A^) 

(2b) 

A^  =  TLAT  -  SLAT 

(la) 

Azimuth  Angle; 

0  =  (ISO/ti)  X  atan2(cp) 

(3a) 

(p  =  cos((t)i)cos({|)2)sin(A7.)/(sin((|)i)  -  sin({|)2)cos(AA,)) 

(3b) 

A?l  =  TLNG  -  SLNG 

(Ib) 

Elevation  Angle: 

v(/  =  (I8O/71)  X  atan(AH/L) 
AH  =  TALT  -  SALT 
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(4) 

(Ic) 


Equations  15  were  given  in  the  introduction  and  show  the  relations  between  modeled, 
observed  and  calculated  angles  that  comprise  geometric  variables. 
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